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A method of monitoring a magnetic field drift of a magnetic resonance imaging apparatus 
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FIELD OF THE INVENTION ^ 

The present invention is related to the field of magnetic resonance imaging 
(MRI), and more particularly without limitation to magnetic field drift compensation. 5>» 

5 BACKGROUND AND PRIOR ART m 
Prior art magnetic resonance imaging apparatus* generate a uniform magnetic q 
field which is also referred to as the main field or the polarising field B 0 . The purpose of the ^ 
magnetic field is that the individual magnetic moments of the spins in the tissue to be 
visualised attempt to align with this polarising field, but precess about it in random order at a 
10 characteristic Larmor frequency which is determined by the gyromagnetic constant of the 
spins and the polarising magnetic field Bo. 

However during operation of a magnetic resonance apparatus the magnetic 
field can vary over time due to various reasons. For example when the magnetic resonance 
apparatus is under heavy duty the passive iron shim pieces of the coils which generate the 
15 magnetic field heat up which leads to a fluctuation of the magnetic field. 

U.S. Patent No. 6,294,913bl shows a method for compensation of variations 
in the polarising magnetic field during magnetic resonance imaging. Monitor signals are 
acquired in an interleaved manner during a scan with the MRI system. Frequency changes 
caused by variations in the polarising magnetic field B 0 are measured using the monitor 
20 signals, and these measured frequency changes are employed to compensate image data 

acquired during the scan. This compensation is achieved by changing the frequency of the RF 
transmitter and receiver to offset the effects of changes in Bo. A disadvantage of this 
compensation method is that special monitor signals are required for the compensation. 

It is therefore an object of the present invention to provide for an improved 
25 method of monitoring and compensating a magnetic field drift. 



SUMMARY OF THE INVENTION 

The present invention provides for a method of monitoring a magnetic field 
drift of a magnetic resonance imaging apparatus which does not require special monitor 



WO 2005/003802 PCT/TO2004/051102 

2 

signals or monitor sensors. Rather the method of the invention can be performed on the basis 
of the magnetic resonance signals which are acquired during magnetic resonance image data 
acquisition, such as by single shot EPI or by a gradient echo sequence. 

In essence the phases of at least two magnetic resonance signals are acquired 

5 an echo time after the corresponding RF excitations. This corresponds to the central k-space 
line which has frequency encoding but no phase encoding. This central k-space is usually 
designated as k y = 0. The difference of two consecutive phase measurements which are 
acquired at a certain time interval provides the shift of the resonance frequency. This enables 
monitoring of the shift of the resonance frequency and compensation of the magnetic field 

10 drift. 

In accordance with a further preferred embodiment of the invention a single 
shot EPI method is used for the magnetic resonance data acquisition. Single shot EPI has the 
advantage that every data acquisition contains the k y = 0 line such that for each data 
acquisition the required phase information can be obtained. 

15' In accordance with a further preferred embodiment of the invention a gradient 

echo sequence is used for the magnetic resonance image data acquisition. During such a 
gradient echo sequence the k-space is scanned along an arbitrary trajectory which typically 
also contains the k-space line k y = 0. For example a full gradient echo sequence has 256 data 
acquisitions one of which is representative of k y = 0. The data acquisition for k y = 0 enables 

20 to obtain the required phase information once for each complete gradient echo sequence. 

In accordance with a further preferred embodiment of the invention extra 
scans along k y = 0 are performed during a gradient echo sequence in order to obtain multiple 
phases during the complete sequence. This enables to monitor a magnetic field drift which 
occurs during a single gradient echo sequence. 

25 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the following preferred embodiments of the invention will be described in 
greater detail by making reference to the drawings in which: 

Figure 1 shows a block diagram of a magnetic resonance apparatus with 
30 magnetic field drift compensation, 

Figure 2 is illustrative of the phase determination for a single shot EPI, 

Figure 3 is illustrative of a flow chart for monitoring and compensating the 
magnetic field drift in the case of single shot EPI, and 
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Figure 4 is illustrative of a flow chart for magnetic field monitoring and 
compensation if a gradient echo sequence is employed. 

DETAILED DESCRIPTION 

Figure 1 shows a magnetic resonance device 1 which includes a first magnet 
system 2 for generating a steady magnetic field, and also several gradient coils 3 for 
generating additional magnetic fields having a gradient in the X, Y, Z directions. The Z 
direction of the co-ordinate system shown corresponds to the direction of the steady magnetic 
field in the magnet system 2 by convention. The measuring co-ordinate system x, y, z to be 
used can be chosen independently of the X, Y, Z system shown in FIG. 3. The gradient coils 
are fed by a power supply unit 4. An RF transmitter coil 5 serves to generate RF magnetic 
excitation pulses and is connected to an RF transmitter and modulator 6. 

A receiver coil is used to receive the magnetic resonance signal generated by 
the RF field in the object 7 to be examined, for example a human or animal body. This coil 
may be the same coil as the RF transmitter coil 5. Furthermore, the magnet system 2 encloses 
an examination space which is large enough to accommodate a part of the body 7 to be 
examined. The RF coil 5 is arranged around or on the part of the body 7 to be examined in 
this examination space. The RF transmitter coil 5 is connected to a signal amplifier and 
demodulation unit 10 via a transmission/reception circuit 9. 

The control unit 1 1 controls the RF transmitter and modulator 6 and the power 
supply unit 4 so as to generate special pulse sequences which contain RF pulses and 
gradients. The phase and amplitude obtained from the demodulation unit 10 are applied to a 
processing unit 12. The processing unit 12 processes the presented signal values (also 
referred to as k-space) so as to form an image by transformation. This image can be 
visualized, for example by means of a monitor 13. 

Processing unit 12 determines a shift Af of the resonance frequency on the 
basis of the phase information obtained from the demodulation unit and provides the shift Af 
to the control unit 1 1 for compensation of the shift. 

The acquisition of the required phase information is explained in more detail 
in the following by making reference to figure 2. First an RF magnetic excitation pulse 200 
is generated by the transmitter coil (cf. transmitter coil 5 of figure 1). Further a slice selecting 
gradient G z pulse 202, a bipolar switched gradient 204, and a phase encoding gradient 206 
are applied. The application of these gradients 202, 204 and 206 for single shot EPI magnetic 
resonance data acquisition is as such known from the prior art. 
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The single shot EPI provides measurement signal 208. The peak of 
measurement signal 208 is the echo time TE after the RF magnetic excitation pulse 200. The 
phase of measurement signal 208 at time TE with respect to the RF magnetic excitation pulse 
200 is determined for the purpose of monitoring a drift of the magnetic field and a 
5 corresponding drift of the resonance frequency. In order to determine the resonance 
frequency drift, if any, the phase of the measurement signal 208 is determined again at a 
subsequent single shot EPI. The phase difference between the single shot EPIs enables to 
calculate the shift of the resonance frequency. This will be explained in greater detail by 
making reference to figure 3 . 

10 The phase information at time TE can be obtained from the measurement 

signal 208 in the time domain or in the frequency domain which requires to perform a Fourier 
transformation on measurement signal 208. The required phase information corresponds to a 
central line of the k-space which is usually referred to as k y = 0, i.e. a scan with frequency 
encoding but without phase encoding. The fact that there is no phase encoding enables to use 

15 drifts of the phase information obtained from the measurement signal 208 for determination 
of the resonance frequency shift. 

Figure 3 shows a flow chart for monitoring and compensating of a magnetic 
field drift in the case of single shot EPI. In step 300 the index d is set to zero. The shift of the 
resonance frequency A/ is also set to zero in step 300. 

20 In step 302 the first single shot EPI is performed and the measurement of the : 

phase 0 a (d) is performed as explained with respect to figure 2. 

The time interval of TE after the first single shot EPI is performed in step 302 
the second single shot EPI is performed in step 304. Again the phase 0 a is measured. 

In step 306 a shift Af of the resonance frequency is calculated by calculating 

25 the difference between the phases 0 a which have been acquired in steps 302 and 304 and by 
dividing the difference by 2 7tx the time interval between the consecutive single shot EPIs. In 
this case this time interval is the echo time TE. 

In step 308 the absolute value of the shift A/ is compared to a threshold value. 
If the threshold value is surpassed a drift compensation is performed in step 310. Preferably 

30 this is done by adjusting the frequency of the RF magnetic excitation pulses to the shifted 
resonance frequency. 
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If the shift 4f does not surpass the threshold no drift compensation is 
required. In step 312 the index d is incremented and the control goes back to step 304 for 
continuous monitoring the shift &f . 

Figure 4 shows an alternative embodiment which uses a gradient echo 
sequence rather than single shot EPI. Step 400 corresponds to step 300 of figure 3. In step 
402 the magnetic resonance data acquisition is started by starting to scan the k-space along a 
given trajectory. The trajectory may or may not cover the central line of the k space which is 
k y = 0. 

In step 404 the scan of the k-space along the trajectory is interrupted in order 
to perform a data acquisition for the central k-space line with frequency encoding but no 
phase encoding, i.e. k y = 0. This way the phase information 0 a is obtained. 

In step 406 the scan along the k-space trajectory is continued. After a certain 
time interval the scan along the k-space trajectory is interrupted again in order to perform 
another data acquisition for k y = 0 to obtain another phase information 0 a . 

In step 410 the shift of the resonance frequency is calculated by calculating the 
difference between the phases determined in steps 404 and 408 and by dividing the difference 
by2;rTE. 

In step 412 the absolute value of the frequency shift Af is compared to a 
threshold value. If the threshold value is surpassed drift compensation is performed in step 
414. In step 416 the index d is incremented and the control goes back to step 406 in order to 
continue the scan along the k-space trajectory. This procedure continues during the entire 
gradient echo sequence for continuously monitoring the shift of the resonance frequency and 
compensating the drift of the magnetic field if necessary. 

The RF pulses of this modified gradient echo MRI sequence are applied at a 
single frequency (for a single slice); the phase measured at the echo (without encoding 
gradients) reflects only the phase difference accumulated between the RF pulse and the TE in 
case that the NMR resonant frequency (due to B0) is not equal to the excitation frequency 
(excluding susceptibility and chemical shift effects). When the magnet B0 is such that the 
NMR resonant frequency is the same as the excitation frequency of the RF pulse then the 
phase error accumulated between RF pulse and TE will be zero or constant from TR to TR. 

When the B0 is such that it corresponds to a different resonant frequency than 
the one the RF pulse is exciting, then the phase error accumulated between RF pulse and TE 
will be a finite value (proportional to the difference and TE). At the peak of each RF pulse, 
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all magnetization is in phase, irrespective of the BO. Only after this time docs the BO have an 
effect on the magnetization. 

So, each RF pulse behaves like a phase reset (as far as the fO measurement 
process is concerned). The relative phase error (relative to zero at the RF pulse peak) at the 
TE increases as BO moves further from the RF excitation frequency. Other mechanisms 
(chemical shift, susceptibility) can contribute to the phase error measured at the TE. 
Assuming that these other contributors to phase error are constant from TR to TR, the change 
in BO can be determined by calculating the difference in phase errors between data acquired 
form two different TR's. 
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